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i-Q ' We provide new values for the model parameters of the covariant constituent quark model (with 

r ^ . built-in infrared confinement) in the meson sector by a fit to the leptonic decay constants and a 

number of electromagnetic decays. We then evaluate, in a parameter-free way, the form factors of the 
1 B(B 3 ) — > P(V) transitions in the full kinematical region of momentum transfer. As an application 

, of our results we calculate the widths of the nonleptonic B s -decays into D~ Df , D* ~ Dj + D~ D* + 

and + . These modes give the largest contribution to Ar for the B s — B s system. We also 

treat the nonleptonic decay B s — > J/if)(j>- Although this mode is color-suppressed, this decay has 
important implications for the search of possible CP- violating new physics effects in B„ — B 3 mixing. 
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I. INTRODUCTION 



(N 
> 

^sO , The study of heavy flavor physics is important due to the unique possibility of determining the Cabibbo-Kobayashi- 
C*~) ■ Maskawa matrix elements. Such studies also provide insights into the origin of flavor and CP-violation. Moreover, 
one of the main purposes of heavy-flavor experiments is to look for new physics beyond the standard model (see the 
recent review [l|). The subject to study are heavy hadrons containing a b- or a c-quark and their weak decays. Note 
that the t quark decays too quickly to form stable hadrons. Recently, time-dependent measurements of CP violation 
in the B s — B s system have become available. In the wake of these measurements, the decay B s — > J/ip cf> has attracted 
much attention from both theorists and experimentalists (see, for instance, Refs. @, Q and references therein). 

The main idea in the theoretical studies of heavy-flavor decays is to separate short-distance (perturbative) QCD 
dynamics from long-distance (nonperturbative) hadronic effects. One uses the so-called naive factorization approach 
which is based on the weak effective Hamiltonian describing quark and lepton transitions in terms of local operators 
that are multiplied by Wilson coefficients (for a review, see Ref. Q). The Wilson coefficients characterize the short- 
distance dynamics and may be reliably evaluated by perturbative methods. The calculation of the hadronic matrix 
elements of local operators between initial and final states require nonperturbative methods. One needs to know how 
hadrons are constructed from quarks. Technically, any matrix element of a local operator may be expressed in terms 
of a set of scalar functions which are referred to as form factors. The so-called QCD factorization and the soft-collinear 
effective theory yield factorization theorems which allow for a systematic description of a given process in terms of 
products of soft and hard matrix elements (we refer an interested reader to Refs. @, @ and Q)- 

A variety of theoretical approaches have been used to evaluate the hadronic form factors. The least model dependent 
among these is the light-cone sum rule (LCSR) approach (see Refs. [§[ and Q). In the LCRS approach one can access 
the form factors in the large recoil (small momentum transfer) region which are then extrapolated to the near- 
zero recoil region using some mo del- dependent pole-type parametrizations. Grinstein and Pirjol have developed a 
systematic approach to the rare decay B — > K*i + i~ in the low recoil region using the heavy-quark effective theory 
framework [lfj| . The low-recoil approach was late studied in detail in Ref. . 

We mention a few other model approaches for the calculation of the form factors. They are based on the study 
of i) Dyson- Schwinger equations in QCD [l2j]; ii) the constituent quark model using dispersion relations [l3[; iii) a 



relativistic quark model developed by Ebert, Faustov and Galkin [14j; iv) a QCD relativistic potential model [15 1 
(see also |l6[); v) a QCD sum rule analysis 1171 ; and, finally, vi) the covariant constituent quark model developed by 
some of us starting with Refs. [Hj], [l!| and [2fl|. It is worth mentioning that the entire physical range of momentum 
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transfer is accessible in the covariant quark model approach used in Refs. [1 81 - 4201 ] and the present paper and in the 

calculations of Refs. [IM1, El- 

The earlier versions of the covariant constituent quark model (for short: covariant quark model (39J) in Refs. fisl - Eol ] 
did not include the confinement of quarks. The nonconfined version had been applied, among others, to the description 
of B and B c transition form factors using a small set of model parameters. In the covariant quark model, meson 
transitions are described by covariant Feynman diagrams with free constituent quark propagation. The ultraviolet 
behavior of the loop diagrams is tempered by appropriately damped vertex functions. A key role in the consistent 
formulation of the model is played by the so-called compositeness condition [2l|,|22j, a corollary of which guarantees the 
correct charge normalization of charged mesons at zero momentum transfer. Since the propagation of the constituent 
quarks is described by free-particle Green's functions, one will encounter on-shell quark production in the case when 
the mass of the bound state exceeds the sum of the constituent quark masses. Therefore, the applicability of the 
covariant quark model in its original version was limited to the cases where mu < m qi + m q2 . This limitation was 
removed later on in Ref. [23j by effectively introducing infrared confinement through the introduction of a universal 
infrared cutoff parameter in the space of loop integrations. This extended the applicability of the covariant quark 
model to all processes involving heavy and light hadrons. The viability of the improved covariant quark model was 
demonstrated in a number of applications to mesonic transitions in Ref. [23| . Later on, this approach was successfully 
applied to a study of the tetraquark state X(3872) and its strong and radiative decays (see Refs. [HIHH). 

Once the parameters of the covariant quark model have been determined the covariant quark model is a very flexible 
tool that can be used to calculate any heavy-to-heavy, heavy-to-light and light-to-light hadron transition. While 
the more model independent approaches usually have to rely on a heavy quark mass expansion the predictions of 
the covariant quark model hold for general mass configurations which are not accessible to the model-independent 
approaches. On the other hand the predictions of the heavy-quark expansion can be recovered by using static 
propagators for the heavy quarks. 

In this paper we use the improved version of the covariant quark model including infrared confinement to evaluate 
the form factors of the B{B S ) — > P(V) transitions in the full kinematical range of momentum transfer. As an 
application of our results we calculate the widths of several B s nonleptonic decays. These are the modes B s — > 
DjDf,D* ~D+ + DjD* + , and B s — > D*~ D* + which give the largest contribution to Ar for the B s — B s system. 
We also treat the color-suppressed mode B s — > J/ip<p. This decay is important for the search of possible CP-violating 
new physics effects in B s — B s mixing. 

Our paper is structured as follows. In Sec. HU we give a brief sketch of the theoretical framework underlying the 
covariant quark model including a discussion of how infrared confinement is set up in the loop integrations. In Sec. lIIIl 
we discuss in some detail how the model parameters of the covariant quark model are determined through a least- 
squares fit to experimental/theoretical data on leptonic decay constants and to eight fundamental mesonic one- and 
two-photon decays. Once the model parameters of the covariant quark model are fixed, the model can be used to 
obtain parameter-free predictions for any transition process involving light or heavy mesons. In Sec. IIV1 we calculate 
the transition form factors of the B and B s mesons to light pseudoscalar and vector mesons which are needed as 
ingredients for the calculation of the semileptonic, nonleptonic and rare decays of the B and B s mesons. In Sec. fVl 
we make use of the calculated form factors to calculate the nonleptonic decays B s — > D s D s and B s —} J/00 
which have recently attracted some attention as explained above. Finally, in Sec. I VII we summarize our findings. 



II. EFFECTIVE LAGRANGIANS, COMPOSITENESS CONDITION, AND INFRARED 

CONFINEMENT 

In this section, we give a brief description of the theoretical framework underlying the formulation of the covariant 
quark model. We first define a nonlocal meson-quark-quark vertex in terms of an effective Lagrangian. We then 
introduce the compositeness condition and discuss its significance. We then, finally, describe how infrared confinement 
is incorporated into the model. This involves a technical discussion of how the one-loop integrations are done, which 
we briefly describe. 

The coupling of a meson H(qiq~2) to its constituent quarks qi and q\ is described by the effective Lagrangian (26l[27j 

AntHqq(a;) = 9hH(x) JdxiJdx 2 F H (x, xx , x 2 )q~2 (x 2 )Tb qi (xi ) + h.c. (1) 

Th is a Dirac matrix or a string of Dirac matrices which projects onto the spin quantum number of the meson 
field H(x). In the present case, the Dirac structures involved are 75 for the pseeudoscalar meson and 7 M for the 
vector meson. The function Fh is related to the scalar part of the Bethe-Salpeter amplitude and characterizes 
the finite size of the meson. To satisfy translational invariance, the scalar function Fh has to fulfill the relation 



3 



Fh (x + a, x\ + a, x 2 + a) = Fjj (x, Xi, x 2 ) for any four- vector a. A specific form which satisfies translational invariance 
is the form 



F H (x,xi 1 x 2 ) = S(x - w\xx - w 2 x 2 )$ H {(%i - x 2 ) 2 ) 



(2) 



where <$>h is the correlation function of the two constituent quarks with masses m qi , m q2 , and the mass ratios 
Wi = m q J(m qi + m q% ). 

The coupling constant gn in Eq. ([TJ is constrained by the so-called compositeness condition originally proposed in 
Refs. [I]], [22j], and extensively used in Refs. [26|, [53] ■ The compositeness condition requires that the renormalization 
constant of the elementary meson field H (x) is set to zero: 



z ff = i-|4K) = o 



(3) 



where H' H is the derivative of the meson mass operator. 

To clarify the physical meaning of the compositeness condition in Eq. ([3]), we first want to remind the reader that the 

1/2 

renormalization constant Zjj can also interpreted as the matrix clement between the physical and the corresponding 
bare state. The condition Zh = implies that the physical state does not contain the bare state and is appropriately 
described as a bound state. The interaction Lagrangian of Eq. ([1]) and the corresponding free parts of the Lagrangian 
describe both the constituents (quarks) and the physical particles (hadrons) which are viewed as the bound states of 
the quarks. As a result of the interaction, the physical particle is dressed, i.e. its mass and wave function have to be 
renormalized. 

In a more familiar setting the compositeness condition Zh — guarantees the correct charge normalization of a 
charged particle at zero momentum transfer. This can be seen by using an identity relating the derivative of the 
free-quark propagator (with loop momentum k + p) with the electromagnetic 7^ coupling to the same propagator at 
zero momentum transfer. The identity reads 



d 



1 



1 



dp^ m q - 



1 



TV 



(4) 



The contribution of the left-hand-side of Eq.Q is normalized due to the compositeness condition, and, therefore, the 
contribution of the right-hand-side is also normalized. 

The condition Zh = also effectively excludes the constituent degrees of freedom from the space of physical 
states. It thereby guarantees that there is no double counting for the physical observable under consideration. The 
constituents exist only in virtual states. One of the corollaries of the compositeness condition is the absence of a 
direct interaction of the dressed charged particle with the electromagnetic field. Taking into account both the tree- 
level diagram and the diagrams with self-energy insertions into the external legs (i.e. the tree-level diagram times 
Zh — 1) yields a common factor Zh, which is equal to zero. We refer the interested reader to our previous papers 



19, 26, 27] where these points are discussed in more detail. 



In the case of pseudoscalar and vector mesons, the derivative of the meson mass operator appearing in Eq. ([3]) can 
be calculated from the one-loop two-point function given by 



1 



d*k 



2p 2 J A-K 2 i 



7 5 5i(fc + wip)j 5 S 2 (k — w 2 p) 
(— k 2 ) |wi tr 7 5 Si(fc + w\p) p'Si(k + wip)^ 5 S 2 (k — w 2 p) 
— w 2 tr 7 5 5*i (k + w\p)"i S 2 (k — w 2 p) p 1 S 2 (k — w 2 p) 
1 ( Vv!Pv\ 1 - d fd A k 



3 \ 9t ^ p 2 



2p 2 ' P dp a J An 2 i 
1 / PnPv\ J_ f d 4 k ~ 2 

2p 2 J 4ir 2 i v 



3 \ 9tn " p 2 



7 $v(-k 2 ) tr *fSi(k + w lP )YS 2 (k - w 2 p) 
(-fc 2 ) jwi tr 7 Al S'i (k + Wip)$Si{k + wip)^S 2 {k - w 2 p) 
w 2 tr 7 p S'i(fc -I- w\p)Y S 2 {k — w 2 p)$S 2 (k — w 2 p) 



}, (5) 
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where (— fc 2 ) is the Fourier-transform of the vertex function $>h{(%i — X2) 2 ), Si(k) is the free-quark propagator 
given by 

Si{k) = —^—s, (6) 
m qi — ft 

and m qi is the effective constituent quark mass m qi . 

For calculational convenience, we will choose a simple Gaussian form for the vertex function fc 2 ). One has 

^ H (-k 2 )=exp(k 2 /A 2 H ) (7) 

where the parameter Ah characterizes the size of the respective bound state meson H. Since fc 2 turns into — k E in 
Euclidean space, the form (JT]) has the appropriate fall-off behavior in the Euclidean region. We emphasize that any 
choice for is appropriate as long as it falls off sufficiently quickly in the ultraviolet region of Euclidean space to 
render the corresponding Feynman diagrams ultraviolet finite. 

The technical details of how the one-loop integrations such as in Eq.© are done can be found in Ref. [23| . Let us 
mention that we use the Schwinger representation to write the local quark propagators as 

00 

S(k) = (m+ JS) f d/3 e- e . (8) 



11 



The loop momentum now appears in the exponent which allows one to deal very efficiently with tensor loop integrals 
by converting loop momenta into derivatives via the identity 

We have written a FORM 28] program that achieves the necessary commutations of the differential operators in a 
very efficient way. 

After doing the loop integration one obtains 

00 

n = J d n /3F(/3 1 ,.. .,/?„) , (10) 


for a given Feynman diagram n, where F stands for the whole structure of a given diagram. For the mass operators 
of Eq. ([5]) one has three propagators, and, thus, one has three Schwinger parameters f3i (i — 1,2,3). For the transition 
form factors to be discussed later on, one has again three propagators leading again to n — 3. 

Next, we briefly describe how infrared confinement is implemented [23J in the quark loops. First, note that the set 
of Schwinger parameters /3i can be turned into a simplex by introducing an additional t integration via the identity 

1= / dt5(t-J2Pi) (11) 
1=1 

leading to 

11 = Jdtt n - 1 Jd n a6(l-jra^F(ta 1 ,...,ta n ). (12) 

1=1 

There are now altogether n numerical integrations: (n — 1) ^-parameter integrations and the integration over the 
scale parameter t. The very large t region corresponds to the region where the singularities of the diagram with its 
local quark propagators start appearing. However, as described in Ref. 23], if one introduces an infrared cut-off on 
the upper limit of the t integration, all singularities vanish because the integral is now analytic for any value of the 
kinematic variables. We cutoff the upper integration at 1/A 2 and obtain 

1/A 2 1 

IT = Jdtt n - 1 Jd n a6(l-J2a l ^F(ta 1 ,...,ta n ). (13) 

i=1 

By introducing the infrared cutoff one has removed all potential thresholds in the quark-loop diagram, i.e. the quarks 
are never on-shell and are thus effectively confined. We take the cutoff parameter A to be the same in all physical 
processes, i.e. the infrared parameter is universal. The numerical evaluation of the integrals have been done by a 
numerical program written in FORTRAN code. 
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TABLE I: Input values for the leptonic decay constants fu (in MeV) and our least-squares fit values. 





Fit Values 


Other 


Ref. 




This work 


Other 


Ref. 




128.7 


1 Qfl A 
lOU.4 


_|_ 
it 




[29, 30J 




198.5 


198 ± 2 


[29J 


Sk 


156.1 


156.1 


± 


0.8 


[29, 30J 




228.2 


227 ± 2 


[29] 


Id 


205.9 


206.7 


± 


8.9 


[29, 30J 




415.0 


415 ± 7 


[29J 


fn s 


257.5 


257.5 


± 


6.1 


[29,30] 




213.7 


217 ± 7 


[29] 


fs 


191.1 


192.8 


± 


9.9 


[2U 


fn* 


243.3 


245 ± 20 


[33] 


fa. 


234.9 


238.8 


± 


9.5 


[31] 


fn; 


272.0 


272 ± 26 


[33] 


fB c 


489.0 


489 


± 


5 


[32] 


h* 


196.0 


196 ± 44 


[33] 


Ip 


221.1 


221 


± 


1 


[29] 




229.0 


229 ± 46 


[33] 



III. MODEL PARAMETERS 



Let us first enumerate the number of model parameters in the covariant quark model. For a given meson Hi these are 



the coupling parameter gn i , the size parameter An i , two of the four effective constituent quark masses 



"i 



(TOi, 



m^, m s , 77i c , m;,), and the universal confinement parameter A. For tin mesons one therefore has 2uh + 5 model 
parameters. The compositeness condition provides nu constraints on the model parameters, which we symbolically 
write as 



fHi (gHi , k Hi ,m qi ,X) = 1 



(14) 



The constraint (1141) can be used, e.g., to eliminate the coupling parameter g B from the set of parameters. The 
remaining parameters are determined by a fit to experimental data. An obvious choice is to fit the model parameters 
to the experimental values of the leptonic decay constants. In the covariant quark model the relevant expressions for 
the pseudoscalar and vector mesons are given by 



N c g P 
N c g v 



d 4 k 
(2tt)H 
d 4 k 



$ P (-fc 2 )tr 



0£S 1 {k + w 1 p) 1 b S 2 {k- w 2 p) 



(27T) 



$ y (-r)tr 



O^Si(k + wip) f< v S2(k - w 2 p) 



= rrivfye 11 



P 



V' 



i P . 



(15) 



where N c — 3 is the number of colors. As before we have Wi — m q J(m qi + m q2 ). Further O p ~ 7 M and Oy = 7^75 . 



The compositeness conditions (1141) and the fit to the leptonic decay constants (|T5j) provide constraint equations 
for 2uh + 5 model parameters. As further constraints on the parameter space, we have decided to fit the model to 
the eight fundamental electromagnetic decays listed in Table HU We refer the reader to Ref. [23j for the details of 
the one-loop calculation of the electromagnetic decays. The results of the (overconstrained) least-squares fit to the 
leptonic decay constants and the electromagnetic decay widths and the corresponding experimental/theoretical input 
values can be found in Tables Q] and Table [Til respectively. The agreement between the fit values and input values is 
quite satisfactory. 

The results of the fit for the values of quark masses m qi , the infrared cutoff parameter A and the size parameters 
A^are given in Eqs. (fTBj) . (fTT)) and (|T5| . respectively. 



7716 



0.235 0.424 2.16 5.09 0.181 GeV 



(16) 



K A K Ad Ad s A b A Bs A Bc A p 



0.87 1.04 1.47 1.57 1.88 1.95 2.42 0.61 GeV 



(17) 
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TABLE II: Input values for some basic electromagnetic decay widths and our least-squares fit values (in keV). 



Process 


Fit Values 


Data [29] 


7T — s> 77 


5.06 x 10" 3 


(7.7 ±0.4) x 10" 3 


Vc -> 77 


1.61 


1.8 ± 0.8 


P — > 7T ± 7 


76.0 


67 ± 7 


U — > 7T°7 


672 


703 ± 25 


K *± 


55.1 


50 ± 5 


K*° -»■ 


116 


116 ± 10 




1.22 


1.5 ± 0.5 


J/0 »7c7 


1.43 


1.58 ± 0.37 



Au A Am A k . A d * A d , A b » A b * 

™ : : (is) 

0.47 0.88 1.48 0.72 1.16 1.17 1.72 1.71 GeV 

The constituent quark masses and the values of the size parameter fall into the expected range. The size parameters 
show the expected general pattern in that the geometrical sizes of the mesons oc A^ 1 shrink as their masses increase. 

The present numerical least-squares fit and the values for the model parameters supersede the results of a similar- 
analysis given in [23| which used a different set of electromegnetic decays in the fit. In the present fit we have also 
updated some of the theoretical/experimental input values. 



IV. TRANSITION FORM FACTORS 



Given the fact that all model parameters have been fixed, the covariant quark model can now be utilized to calculate 
any given decay process in a parameter-free way. As a first application we calculate the form factors describing the 
transitions of heavy B{B S ) mesons into light mesons, e.g. B,B S — > ir, K, p, K*, <fi. These quantities are of great 
interest due to their applications in semileptonic, nonleptonic and rare decays of the B and B s mesons. They have 
been calculated within the LCSR approach in the region of large recoil (small momentum transfer) and have been 
extrapolated to the low recoil region. Our approach allows one to evaluate the form factors in the full kinematical 
range including the near-zero recoil region. 

Below, we list the definitions of the dimensionless invariant transition form factors together with the covariant 
quark model expressions that allow one to calculate them. We closely follow the notation used in our papers [l9[ . 



(P[, 3q2 ](P^\^O^ qi \P lq3gi] ( Pl )) 
N c gp9p' I — — ®i 



{2tt)H 
F + (q 2 )PV + F-(q 2 )q>* 



( - (fc + wis) 2 ) $p> ( - (* + ™ 23 ) 2 )tr 



0»S 1 (k + p 1 ) 1 5 S 3 (k) 7 5 S 2 (k + P2 ) 



(19) 



(%. }q2 ](P2)\q2(cr^q,)qi\P [q3qi] (pi)) 
, d 4 k 

N c gp g P 
i 



.,_ , H $p ( - (* + ^i 3 ) 2 ) §p> ( - (k + w 23 ) 2 )tY 



a » v q v Si (k+ Pl ) 7 5 S 3 (k) 7 5 S 2 (k + p 2 ) 



mi + m 2 



fP^-q-Pq") F T (q 2 



(20) 
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(V(p 2 , e2)[g 3 g 2 ] I 92 qi \ P[ q3qi ](pi)) 
= N c gpg v [jll^$p(-(k + w 13 ) 2 ^ $ v (- (fc + w 23 ) 2 )tr 



O^i(fc+Pi) 7 5 S 3 (£0 ^(fc+pa) 



mi + 771 2 



(2tt) 4 7 

(-<T P • g^(9 2 ) + P M A + (q 2 ) + q»P v A_{q 2 ) + 7 ^ P a q fj V(q 2 )) , (21) 



(V(p 2 , e 2 ) te92 ] I 92 (<t mi/ 9v(1 + 7 5 )) 9i I % 9l ](Pl)) 

N c gp gv J $p ( - (fc + wis) 2 ) $v ( - (fc + w 23 ) 2 ) tr 



(<r ""?„(! + 7 5 )) Sx(k +Pi) 7 5 5 3 (*) ^ S 2 (k + 



P2) 



= 4 {-(g» u - q»q" /q 2 )P ■ qa a {q 2 ) + {P» P» - q» P v P ■ q/q 2 )a+{q 2 ) +ie^ P a qpg{q 2 )) . (22) 

We use P = pi + p-2 and q = p\ — p 2 and the on-shell conditions e\ • P2 = 0, pf = mf. Since there are three quark 
species involved in the transition, we have introduced a two-subscript notation Wij — m qj /(m qi + m qj ) = 1, 2, 3) 
such that Wij + Wji — 1. The form factors defined in Eq. (|22l) satisfy the physical requirement cto(0) = a+(0), which 
ensures that no kinematic singularity appears in the matrix element at q 2 — 0. For reference, it is useful to relate the 
above form factors to those used in Ref. Q. The relations read 

*+ = /+, P- = -^Sr^(/+-/o), F T = f T , 

9 

A = r2l±^ Al , A + = A 2 , A _ = 2m 2 (m 1 + m 2 ) y = y 

777l — ffl2 q 

9 2 

ao = P2 , g = T t , a + = T 2 H 5 2 T 3 • (23) 

— 777-2 

The form factors (|23p satisfy the constraints 

4,(0) - A 3 (0) (24) 

2777 2 A 3 (9 2 ) = (7771 + 777 2 )^l(g 2 )- (7771 - 777 2 )A 2 (g 2 ). 

In Figs. QUI we plot our calculated form factors in the entire kinematical range < q 2 < q^ax- For comparison, 
we also show the results obtained from the light-cone sum rules analysis [1| . In Table IIII| we collect our predictions 
for the form factors at the maximum recoil point q 2 — and provide a comparison with results obtained within other 
approaches. The figures and tables highlight the wide range of phenomena accessible within our approach. 

As was suggested in Ref. [ll|, one can check how well the quark-model form factors satisfy the three low recoil 
relations derived in Ref. [l(| involving the pairs of form factors (Ti, V), (T 2 , A{), and (T 3 , A 2 ). In Fig. [SJ we plot the 
ratios 

r Tliq2) r T2{q2) r q2 T3iq2) m) 

which in the heavy-quark symmetry limit and at low recoil, should all be of the order 1 — (2a s /(37r) In (/j/mj), i.e. 
close to one. Figure [5] shows that, similar to the extrapolated LCSR form factors, the covariant quark model form 
factors satisfy the low- recoil heavy-quark symmetry relations reasonably well for Pi and P 2 but not for P 3 . Note 
that the q 2 scale has changed in Fig. [5] 

It is interesting to compare the q 2 behavior of the B — ir transition form factors calculated from the three-point 
one-loop diagram with the q 2 behavior of the vector-dominance model (VDM). For example, in a monopole ansatz 
for the form factor F^ 7T {q 2 ) one would have the VDM q 2 behavior 



p Btt (J2\ _ ^+ 7r (°) 
^VDMW I ~ —2 ~2 ' 



m B , — q^ 

where the pole mass is given by the mass of the lowest- lying vector meson state B* . The two curves are plotted in 
Fig. |6l One observes a strong rise of the VDM form factor towards the larger q 2 values close to the position of the 
B* pole. A similar rise is observed for the quark-model form factor. It is quite intriguing and gratifying that the 
quark-model form factor is able to emulate the pole-type behavior of the VDM form factor including even the correct 
scale m B of the pole-type enhancement. 
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TABLE III: q 2 — results for the transition form factors in various model approaches. 



This work LCSR-1 [8] LCSR-2 [9] DSE [12] QCD SR [17] DQM [13] RQM [14] RCQM [18] 



Ff"(0) 


0. 


.29 


0.258±0.031 


0.25±0.05 


0. 


,24±0 


.05 


0.24±0.03 





.29 


0.22 


Ff K (0) 





.42 


0.335±0.042 


0.31±0.04 





,30±0 


.06 


0.254=0.03 





.36 




7^(0) 


0. 


.27 


0.253±0.028 


0.21±0.04 


0. 


,25±0. 


.05 







.28 




F^ K (0) 


0. 


.40 


0.359±0.038 


0.274=0.04 


0. 


32±0 


.06 


0.144=0.03 





.35 




V Sp (0) 


0. 


.28 


0.324±0.029 


0.32±0.10 


0. 


.314=0 


.06 







.31 


0.30 


1/ Sif (0) 


0. 


.36 


0.412±0.045 


0.39±0.11 


0. 


.374=0 


.07 


0.474=0.03 





.44 






0. 


.32 


0.434±0.035 


















Af p (0) 


0. 


.26 


0.240±0.024 


0.24±0.08 


0. 


.244=0. 


.05 







.26 


0.27 


Af K * (0) 


0. 


.33 


0.290±0.036 


0.30±0.08 


0. 


.294=0. 


.06 


0.374=0.03 





.36 




Af s ^(0) 


0. 


.29 


0.311±0.029 


















Af"(0) 


0. 


.24 


0.221±0.023 


0.21±0.09 


0. 


.25±0 


.05 







.24 


0.28 


A^*(0) 


0. 


.32 


0.258±0.035 


0.26±0.08 


0. 


30±0 


.06 


0.404=0.03 





.32 




A***(0) 


0. 


.28 


0.234±0.028 


















if (0) 


0. 


.25 


0.268±0.021 


0.28±0.09 


0. 


26±0 


.05 







.27 




if*" (o) 


0. 


.33 


0.332±0.037 


0.33±0.10 


0. 


30±0 


.06 


0.194=0.03 





.39 




if^(o) 


0. 


.28 


0.349±0.033 



















V. NONLEPTONIC B s DECAYS 

As a second application we consider the two-body nonleptonic B s decays B s — > D s D s and B s — > J/0</>, 
which have recently attracted some interest. The modes DjDf, D*~D+ + Dj D* + , and D*~D* + give the largest 
contribution to Ar for the B s — B s system. The mode J/0</> is color-suppressed but it is interesting for the search of 
possible CP-violating new physics effects in B s — B s mixing. 

It is convenient to express all physical observables in terms of helicity form factors H m . This will result in very com- 
pact rate expressions. Furthermore, in the case of the two-vector meson decays P — > VV, the helicity representation 
is quite convenient since one can then easily calculate the helicity composition of the rate T^, F_, T + . 

The helicity form factors H m can be expressed in terms of the invariant form factors of Ref. H in the following way 




(a) Spin S = 0: 



Hi 



{{m\ - ml) F+ + q 2 F_} 



(26) 



Hq = 2mijp a | ^ 



(b) Spin 5= 1: 
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FIG. 1: Our results for the form factors appearing in Eqs. (|19|) and (I20|l - Left panel, B — ir— transition; and right panel; 
B — A'transition. For comparison, we plot the corresponding LCSR curves from Ref. 



H t = ) milp ^ {( m 2 - m 2 )(A + -A ) + q 2 A_}, 
mi + m 2 m 2 y/q 2 



H± = ; {-(m 2 1 -m 2 2 )A ±2m 1 \p 2 \V} 1 (27) 

mi + 7712 

H = \ ^- T ={-{m 2 1 -ml)(m 2 1 -m 2 2 -q 2 )Ao + 4mj\p 2 \ 2 A + }, 

mi + 77i2 2m 2 V9 " 

where |p>2 1 = X 1 ^ 2 (m 2 , m 2 , q 2 ) / (2 mi) is the momentum of the outgoing particles in the rest frame of the decaying 
particle. 

The effective Hamiltonian describing the B s nonleptonic decays is given by (see Ref. [1]) 

v 2 i=i 

Ql = {Ca 1 ba 2 )v-A(Sa 2 Ca 1 )v-A, Q2 — (c ai 6 0l )v~A, («a 2 Ca 2 )v~A, 

Q3 = {Sa- L ba 1 )v-A(Ca 2 Ca 2 )v-A, Qi — {s ai b a2 )v -A{c a2 C ai )v -A, 

Qb = (Sa 1 ba 1 )v-A(c a2 Ca 2 )v+A, Qq = {s ai b a2 )v—A{c a2 C ai )v+Ai (28) 
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10 
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FIG. 2: Our results for the form factors appearing in Eqs. (|21|) and (122 p for B — p transition. For comparison, we plot the 
corresponding LCSR curves from Ref. ||. 



where the subscript V — A refers to the usual left-chiral current = 7 M (1 — 7 5 ) and V + A to the usual right-chiral 



one O't — 7 M (1 + 7 5 ). The a, denote color indices. We calculate the nonleptonic _B s -decay widths by using naive 



factorization. In this paper we consider the following nonleptonic decays of the i? s -meson: 

B s (p)^D;( qi )Df(q 2 ), Dj(qi)D* s + (q 2 ,e 2 ), D* s ~ ( qi , e 1 )Df(q 2 ), D* g - (q 1 ,e 1 )D* + (q 2 ,e 2 ), and Jft( qi , ex) cj>(q 2 , e 2 ). 

The widths can be conveniently expressed in terms of the helicity form factors and leptonic decay constants. In the 
case of the color- allowed decays B s — > D\ ' D s one has 



r(s s 

T(B S - 
T(B S - 

r(B s ^, 



Here, \ ( c s) 



D-Dt 



d-d: 



^D* s -Df) 
D*-D*+) 



off fPS p B s D B ( 2 \ 



rj en f 

JD S m S 



F 1^2 I r\(s)l2 ( s-,e& m f uB s D s(rn 2 

7T m B \ 

Gf |q 2 | r,( s ) l2 ( r eff m f ttB s D* s( 2 \ , o^cff ragjq 2 | fPS ^B s D' bI _2 \ 

2 



Gf 


|q 2 | 


167r m| 


Gf 


|q 2 | 


16tt 




Gf 


|q 2 | 


16n trig 


Gf 


|q 2 | 


16tt 





F |q 2 | rx ( B )!2 / ^eff J* Ff B s D s( 2 \ 

-^-l A c I ^ C 2 m D*jD*B Q {m D .) 



i=0,± 



(29) 



x/2 



# |Vrf,V^|. The Wilson coefficients appear in the combinations C 2 eff = C2 + £ Ci + C4 + £ C3 and 



C 6 = Cq + £65. where terms multiplied by the color factor £ = 1/N C will be dropped in the numerical calculations 
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B-K*: A,(q 2 ) B-K*: V(q 2 ) 



0.6 




q 2 (GeV 2 ) q 2 (GeV 2 ) 



FIG. 3: Our results for the form factors appearing in Eqs. (|21[) and (|22|l for B — K* transition. For comparison, we plot the 
corresponding LCSR curves from Ref. 0. 



according to the 1/N C — expansion. The annihilation channels that also contribute to the above color-allowed decays 
will be neglected since they are color and form factor suppressed. 
The width of the color-suppressed B s — > J/tp <f> decay is written as 

T(B S -+ J/**) = ^^m 2 (Cf + Cf) 2 {m m f m f £ (H?^{m* m ))\ (30) 

B ° i=0,± 

where we have combined the Wilson coefficients into C£ s = d + £ C 2 + C 3 + £ C 4 and C 5 eff = C 5 +£,C 6 .. 
For the Cabibbo-Kobayashi-Maskawa-matrix elements, we use the values from Ref. [29j: 

\V ud \ \V US \ \V ub \ \Ved\ \Vcs\ \Vcb\ 

(31) 

0.974 0.225 0.00389 0.230 0.975 0.0406 

For the values of the Wilson coefficients we take [34| 



C/~< s~i f-1 S~1 

1 <^2 L-3 W ^5 <^6 

-0.257 1.009 -0.005 -0.078 0.000 0.001 



(32) 
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FIG. 4: Our results for the form factors appearing in Eqs. (|21|) and (|22[) for B s — tj> transition. For comparison, we plot the 
corresponding LCSR curves from Ref. 



evaluated to next-to-next-to leading logarithmic accuracy in the MS (NDR) renormalization scheme at the scale 
H = 4.8 GeV 

We also need the values of the B s — <fi transition form factors evaluated at q 2 = m 2 ]/^, ■ They are given in Table IIVI 
where we compare our results with corresponding results of Ref. [3[. The agreement for the form factors Ai(mj,^) and 
^(m^) is satisfactory. Our value for the form factor V(m 2 j^) is somewhat smaller than the one found in Ref. [3j. 

TABLE IV: The relevant B s — <j> form factors at q 2 — m 2 j^ calculated in this work. For comparison, we give the results of 
Ref. 0]. 





This work 


Ref. [3] 




0.37 


0.42±0.06 


Ai{m 2 m ) 


0.48 


0.38±0.06 




0.56 


0.82±0.I2 



In Table fVl we give our results for the branching ratios. One can see that there is good agreement with the available 
experimental data. 

We finally give our results on the helicity fractions in the two decays B s — > D* ~D* a + and B s — > J/tp4>. The helicity 
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FIG. 5: Our results for the ratios of the form factors appearing in Eq. ((25} for B — .ft"* —transition. 

B-Jt form factor from loop and B -monopole 




q (QeV) 



FIG. 6: Comparison of the B — n— form factors obtained from the covariant quark model and from a VDM-monopole ansatz. 



fractions of the nonleptonic B s 
Wo\ 2 



VV rates are defined as 



\H±\ 



\Ho\ 



H 4 



\Ho\ 



H, 



2 Hq 



(33) 



Note that we have normalized the partial helicity rates to the total rate such that one has (f l + f _ + f +) = 1- 
For B s -> D*-D*+ we find (f L ,f_,f+) =(0.549,0.366,0.0847) and for B s ->■ Jfip<f> we find (0.420,0.552,0.0272). 
The hierarchy of partial helicity rates > T_ > T + seen in the decay B s — > D*~D* + is expected for tree-level- 
dominated nonleptonic decays using simple on-shell quark model arguments. One finds that, at the leading order 



of mi = m,B s , the partial rate T_ is helicity-suppressed by the factor 4q 2 /mf with q 2 



and the partial rate 



r+ is further chirality suppressed by the factor rri^/m\ with mi 



in addition to the helicity suppression |3j 



|38| |. Using the qualitative suppression factors one finds (0.583,0.361,0.056) for the helicity fractions in the decay 
B s — > D*~D* + which is remarkably close to the results of the full calculation. For the process B s — > Jfycj) with a 
larger g 2 -value of q 2 = m 2 j^ the helicity suppression is no longer in effect since now Aq 2 jm\ = 1.332. One now obtains 
(0.420, 0.560, 0.020) for the helicity fractions which again is remarkably close to the results of the full calculation. One 
has an inversion of the hierarchy for the longitudinal and transverse-minus rates for B s — > J/i/j4> m as much as one 
now has 1^ < T_. Experimental numbers on the partial helicity rates exist only for the decay B s — ¥ J/iJj4> given by 
f l = 0.541 ± 0.017 and f j_ = 0.241 ± 0.023 29]. Our calculated longitudinal rate can be seen to be off by several 
standard deviations. In order to be able to compare with the experimental transverse rate T±_ one needs to use 
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TABLE V: Branching ratios (%) of the B s nonleptonic decays calculated in our approach. 



Process This work Exp. data [29] 

B s ->■ D~D+ 1.65 

B s D-D*+ + D*-Df 2.40 

B S ^D*-D*+ 3.18 

B a ->■ J/00 0.16 



r ± oc \A±_\ 2 = \H+ - iJ_| 2 /2. For B s -> D*~D*+ we find t± = 0.0493 and for B s -> J/ip<p we predict t± = 0.167. 
Again we are off the experimental result by several standard deviations. 

VI. SUMMARY 

We have given a brief sketch of the theoretical framework underlying the covariant quark model, including a 
discussion of how infrared confinement is incorporated in the model. We have discussed in some detail how the 
model parameters of the covariant quark model are determined through a least-squares fit to experimental/theoretical 
data on the leptonic decay constants and eight fundamental mesonic one- and two-photon decays. Once the model 
parameters of the covariant quark model are fixed the model can be used to obtain parameter-free predictions for any 
transition process involving light or heavy mesons. 

In the present paper, we have calculated the transition form factors of the heavy B and B s mesons to light 
pseudoscalar and vector mesons, which are needed as ingredients for the calculation of the semileptonic, nonleptonic, 
and rare decays of the B and B a mesons. Our form factor results hold in the full kinematical range of momentum 
transfer. We have provided a detailed discussion of how the covariant-quark-model form factors compare with the 
corresponding LCSR form factors. 

We have finally made use of the calculated form factors to calculate the nonleptonic decays B s — > D s ' D s ' and 
B s — > J/00, which have been widely discussed recently in the context of B s — £? s -mixing and CP violation. We have 
also presented results on the helicity composition for the decays B s — > VV. Further application of our form factor 

results are envisaged, such as the calculation of the penguin-dominated decay B s — > K ( *> K 
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